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Abstract 

The quantum dynamics of the Josephson junction system in the computa- 
tional subspace is investigated. A scheme for the controlled not operation 
is given for two capasitively coupled SQUIDs. In this system, there is no 
systematic error for the two qubit operation. For the inductively coupled 
SQUIDs, the effective Hamiltonian causes systematic errors in the computa- 
tional subspace for the two qubit operation. Using the purterbation theory, 
we construct a more precise effective Hamiltonian. This new effective Hamil- 
tonian reduces the systematic error to the level much lower than the threshold 
of the fault resilent quantum computation. 

I. INTRODUCTION 

It has been shown that, a quantum computer, if available, can perform certain tasks 
much more emciently than a classical Turing machine |ï]-f|]. The realization of the bàsic 
constitute of quantum computer, fault tolerate quantum lògic gate, is a certral issue in the 
subject. One can make a fault tolerate quantum computation through the quantum error 
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correction. It has been pointed out that, if the error rate of each operation is smaller than 
certain threshold, the error correction method works for arbitrary large scale computation. 
The threshold is estimated to be 10~ 6 . 

Recently, it is proposed [|]-[7|] to implement the the quantum gate by súper conducting 
interference device(SQUID). The single qubit operation has been experimentally demon- 
strated in ref @]. This implement ation has a unumber of advantages. One can address the 
single qubit ( Josephson Junction) instead of the bulk material in NMR. One can even con- 
trol the interaction between the two qubit system through the external parameters( voltages 
and magnètic flux). However, strictly speaking, the Josephson junction is not a two level 
system. Only when we choose certain specific parameters we can approximately have a two 
level computational subspace. In general there is a small transition probability between the 
computational spapce and the outside space, this causes the leakage ||. The leakage error 
can be minimized by the specifically designed farbrication ||. Also, one may detect and 
correct the leakage error by taking the meassurement on the space outside computational 
space. Besides the leakage, the approximate effective Hamiltonian causes another type of 
error, the phase shift error. If the rate of all accumulatable errors are lower than certain 
threshold |T7|], a large reliable quantum computation can be done, in principle. Because in 
such cases we can always take certain appropriate error correction procedure to reduce the 
error exponentially ||17|| . In this paper we analyse the phase shift error caused by the approx- 
imate effective Hamiltonian. We give a new effective Hamiltonian through the purtabation 
method. Using this new effective Hamiltonian, the phase shift error rate is lower than the 



quantum computation threshold fï7| . 

Consider a superconducting electron box formed by a symmetric SQUID( see fig. HJ), 
pieced by a magnètic flux <3> and with an applied gate voltage V x . The device is operated in 
the charging regime, i.e., the Josephson couplings Ej are much smaller than the charging 
energy E c h- Also a temperature much lower than the Josephson coupling is assumed. The 
Hamiltonian for this system is |T3]|||7| 



2 



H = E ch (n - n x ) 2 - Ej{$) cos(0), (1) 

Ej(<È) = 2E J0 cos í 71 "^-) j -^c/i is the charging energy and n x can be tuned by the applied 
voltage V x through V x = 2en x /C x ( see figure [![). The phase difference across the junction 
X and the cooper pair number n canonical conjugate variables [6,n] = i. $o = h/2e is the 
quantum of flux. So here n x and $ can be tuned externally. 

If the value of parameter n x is close to 1/2, the energy gap between the ground state |0) 
and the first excited state |1) is much smaller than the gaps among any other states. Thus 
the basis |0), |1) approximately make a computational subspace. The transition between the 
states in the subspace and the state outside the subspace is small. Through calculation of 
the matrix element i.e., (niH\n2){ ni,n 2 = 0, 1), the matrix form of the Hamiltonian H in 



the subspace |0), |1) is 



H e = -~[B x a x + B y a y + B z a z ), (2) 



where (B x , B y , B z ) = (Ej(<&), —Ej(&), E ch {\ — 2n x ) and cr x ^^ are Pauli matrices. Intutively, 
we may regard this H e as the effective Hamiltonian in the computation subspace. We can 
take the rotating operation throhgh the time ecolution property of this effective Hamiltonian 

However, this H e may not be the best choice in simulating the real time evolution. To 
reduce the systematic error, one can use another form of the effective Hamiltonian which 
simulates the time evolution in the computational subspace more precisely. In this paper, 
we calculate the matrix element e~ lHt by the purterbation expansion. We know the time 
evolution operator generated by H can be expressed in the form e~ lHt = J2n e~ tE ' nt \^ ' n ) (\P r n |, 
we still have difficulty in obtaining the precise time evolution properties for state |0) and |1). 
Here E' n and \^ n ) are eigenvalue and eigenstate of H, respectively. \^ n ) can be expressed in 
Mathieu function, which is the infinite summation of trignometric functions. But in the real 
calculation, we have to use the truncated Mathieu function. As we will show it latter, after 
taking the first non-zero perturbation term into consideration, the error is much smaller 



than the threshold for fault resilent quantum computation [0. In the following sections, 



after the we study the single qubit case by the perpurbation method, we will investigate the 
two qubit dynamics, which is the heart of the elementary quantum lògic device. We will 
study both the capasitively coupled SQUIDs and the inductively coupled SQUIDs. For the 
former one, we find the effective Hamiltonian in the subspace works exactly. For the latter 
one, we give a new form for the effective Hamiltonian so that it works more precisely in the 
computational subspace, especially, the error rate is much lower than the threshold for a 
large, reliable quantum computation. 

II. TIME EVOLUTION FROM THE PERTURBATION TREATMENT FOR ONE 

QUBIT 

We first write the Hamiltonian H in the following equivalent form 



H = E ch J2(n - n x ) 2 \n)(n\ - £(|n)(n + 1 + \n + l)(n\). (3) 

1 n=0 



Furthermore, we decompose the Hamiltonian into two parts as 

H = H + H c 

and 



H = E ch ^(n - n x ) 2 \n)(n\ - ^(|0)(1| + |1)(0|) = H e + E ch ]T(n - n x ) 2 \n)(n\, (4) 

n * n >l 



E 00 

^e = ^E(H(^+ 1 + l^+ 1 )H)· (5) 

We regard this H c as the perbation term. H part can be solved exactly. The first two 
eigenvalues(_Eo°' ) and e[ ^) and first two eigenstates(|0 o ) an d \4>i)) are just that of H e . The 
rest eigenvalues and eigenstates are are just that of operator E ch J2 n ( n ~~ n x) 2 \n)(n\. They 
have already been denoted as E n and \n), respectively. If we ignore H c , this H Q is identitical 
to H e in the subspace of |0), The perturbation mordification comes from H c . We do 
the perturbation calculation to the second order of for the eigenvalues and to the first order 



to eigenstates here. (For more exact result, one can do it similarly through higher order 
approximation.) We heve the following results: 



E ch (l-2n x + 2n 2 x )-^EÏ h (l -2n x )+Ej 
E = 7T H A ^o, (6) 



E ch {l-2n x + 2nl) + ^El h {l-2n x )+E] 
E\ = h AL·i {() 

E 2 E 2 

Here AE X = . p „ — r and AE = -r^r 1 — Take the perturbation result above into consid- 
eration, we obtain the time evolution operator in the computational subspace 

^) = Ee-^|0n)(0n|. (8) 
ra=0 

Therefore, the new effective Hamiltonian with this timeevolution property is 

H' e = £ o |0qX0o| + E^i^l =H e + AE \<f) )(<p \ + AE^fa) (fa] (9) 

This new effective Hamiltonian gives the evolution operator as Eq. (^j). If n x is close to 
1/2, the modification in the effective Hamiltonian is insigficant. Especially, if n x = 1/2, 
H' e is same with H e , up to a constant term. But they may differ obviously when n x is far 
from 1/2. With this new effective Hamiltonian, if we omit the leakage error, the distance 
between the state under real time evolution (\ip r (t))) and the state given by evolution of Eq. 

hi E 4 

(i) i s smaller than Typically, if Ej/E c h = 0.02, the error is is between 10 

ch 

to 10 -7 , smaller than the threshold for resilent quantum computation, 10~ 6 to 10~ 5 HÏ7 



However, if we use the old effective Hamiltonian which is the first order perturbation in 
the computational subspace, the error rate is in the magnitude order of 10~ 4 , larger than 
the threshold. However, the error rate can be larger than this threshold if we use the old 
effective Hamiltonian H e , provided that n x is far from 1/2. If the value Ej/E c h is not so 
small in certain case, we can take a higher order purterbation calculation. 



III. TWO QUBIT DYNAMICS 
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A. capasitive coupling 

For the capasitively coupled SQUID, the interacting Hamiltonian is 

Hi = S(n 1 -n Xt i)(n2-n Xí2 ), (10) 

5 = 2E c f l Ex / C. Now we regard qubit 1 as the comtroll bit and qubit 2 as the target bit. Our 
goal is to make a conditional flip to qubit 2, i.e., if state of qubit 1 is |0), nothing happens 
to qubit 2; if state of qubit 1 is |1), then qubit 2 is flipped. To make the C-NOT gate, we 
set n Xt i to and $i = $o/2 for qubit 1. For qubit 2 we set the magnètic flux $ 2 = so 
that Ej = 0; and n Xi 2 = 1/2. Just wait for a period of 

Aí = n/5. (11) 
We obtain the following conditional unitary transformation for qubit 2 

U 

if state of qubit 1 is |0); and 

U 

if state of qubit 1 is |1). This is equivalent to C-NOT gate through Hardmard transformation. 
Since Ej is set to here, so the above two qubits operation is exact. Thus if we use the 
capasitively coupled SQUIDS, there is no phase shift or leakage error caused by the the states 
outside the computational subspace in the two qubits operation. 

B. Inductively coupled case 

The C-NOT gate of inductively coupled SQUID system is given in by using the 
effective Hamiltonian in the computational system. Here we use the time evolution operator 



(i \ 



(12) 



(13) 
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in the whole system and then to calculate the time t at which the fidelity between time 
evolution operator and a C-NOT matrix has the maximum value. 
The Hamiltonian for the two interacting SQUIDs is 



H = H 1 + H 2 - -LíEjt^) sinGi + £ J2 ($ 2 ) sin0 2 ) 2 

We know sinOj = (\ni)(rii + 1| — \n>i + l)(n|)/(2i). We denote 

Ej\Ej2 



(14) 



H = H ®l + l®H + 



2Er 



[(|0)(1|-|1)(0|)®(|0)(1|-|1)(0|)] 



(15) 



and 



H r II l I ll- Kj ^' n I E|n)(n + l| 

2L· l \ n ^ 



// -+- l)(n| Cg> \n)(n + 1| — \n + 



Ej2 /+oo \ cj2 

+ ^(^|n)(n + 2| + ]n + 2)(n| ®1 + 72 



' +oo 



^|n)(n + 2| + |n + 2)( 



n\ 



(16) 



^'z- V-oc / E L 

Previously, the Hq above was used for the effective Hamiltonian in the computational sub- 
space for two qubits operation ^-0] It is easy to see the total Hamiltonian for the two 
SQUIDS is 



H = H + H 



(17) 



We can regard Hp as the perturbation term. In the two qubit operation, we can always set 
n x i = n x2 = 1/2. Thus H e is simplified to Ej($) (|0)(1| + |1)(0|). Under this condition we 
have 



H Q = U 



Ejx + E J2 

Eji - E J2 L·i 

EjiEj 2 j£ 



E. n Ej2 



E,j\Ej2 
\ El 



El 







E.JlEj2 Q 

J2 — Ej\ 

—Ejx — E J2 



\ 



(18) 



where we have used the basis of 1 00 >, 1 1 >, 1 10 >, |11 > and U — e ± ay ®e . Explicitly, 

d-i\ 

We can obtain the eigenvalues and eigenstates of Hq exactly. The 



e ~ Gy 



V2 
2 



1 1 



/ 



i? 2 E 2 i E 2 _E 2 

four eigenvalues are ±J (E J1 + E J2 ) 2 + J E i J2 , ±J (E J1 - E J2 ) 2 + J1 EÍ J2 , corresponding to 



the following four eigenstates respectively 
|^oo) = f/(cos||00) + sin||ll) N 



[ x+ (|oo) + |ii)) + x-(|oi) + |io))], 



and 



Here 



|Voi> = U (cos||01) + s in||10>) = 77+(-|00) + |11)) + 77-flOl) - |10)), 
l^io) = U fsin||01) - cos||10>^) = 77 + (-|00) + |11» - V -(\Q1) - |10>, 



|^i> = C/ (sin||00) - cos ||11>) = 1 [ X+ (|00> + |11» - x _(|01> + |10»] 



X± 



sin|±cosf,77± =sinf icosf ,cos^ = (E^+E^)/^ (E n + E J2 f + (^èf 22 ) ,sin0i = 
£n£ J2 / (el^Ej! + E J2 y + (^èf^) 2 ) , cos fl 2 = - Ej 2 / ] /(Ej l -Ej 2 y + (^ E f^)' 



and sin# 2 = -E n E J2 / J(E. n + E J2 f + ^ 



E.nE, 



E T 



Now we take the first non-zero modification in the perturbation to the 4 states. The modi- 
fications to the 4 energy levels above are 

A + ,0,0A_, 



where 



and 



AE 



ii 



8E, 



ch 



E J1 + E J2 + 



1 ( EjiE J2 

Et 



+ 



E J1 + E J2 
Et 



(19) 



AE, 



oi 



8E ch 



E J1 + E J2 



1 ( Ej\E J2 
Et 



2 ( Ej 1± Ej 2 \ 2 ' 

3 1 3£ L i 



(20) 



Up to an unimportant constant term, the new effective Hamiltonian is 
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H E = H + A + |Voo)<^oo| + A_|Vn)(^n| (21) 

Suppose initially we have a state after time t, the state in the computational subspace 
is p(i) = Tle~' lHt \ty i) i\e lHt Tl, l·l — \ kj){kj\ and k,j can take vàlues of and l. However, 
if we use the effective Hamiltonian, the state we supposed in the computational subsspace 
is p e (t) = T\.e~ lHEt yS> i) i\e lHEt T\.. To estimate the error rate caused by the phase shift error 
due to the effective Hamiltonian, we just calculate distance between p(t) and p e (t). Suppose 
the vàlues of Eji, Ej 2 are close, El is not larger than Ej\ or Ej 2 . After calculation we 
know, if we use the Ho for the effective Hamiltonian, the magnitude order of the distance 

E 2 

can be w ü . However, if we use the modified one, H E for the effective Hamiltonian, the 

ch 

E A 

error rate is reduced to the magnitude order of 10 ~t . This is to say, if Ho is used as the 

ch 

effective Hamiltonian, the error rate is larger than threshold for the fault resilent quantum 
computation, 1CT 6 to 1CT 5 , provided that Ej/E ch is larger than 0.01. However, if we use 
H E , the one proposed in this paper, the error rate is much smaller than the threshold value 
for the fault resilent quantum computation. 

In summary, the properties of the effective Hamiltonian of the Josephson junction system 
in the computational subspace is investigated. For capasitively coupled Josephson junction 
system, there is no systematic error due to the effective Hamiltonian in the computational 
subspace for the two qubits operation. If the inductively coupled SQUIDs are used, the 
effective Hamiltonian in the computational subspace causes phase shift error. But the this 
effective Hamiltonian can be modified to a more exact form. Our new effective Hamiltonian 
given by the perturbation theory reduces the systematic error to a range much lower than 
the threshold of the fault tolerate quantum computation. 
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FIGURES 
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FIG. 1. SQUID with symmetric Josephson junctions 
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FIG. 2. The capasitively coupled SQUIDs. 
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